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ABSTRACT 



H3 has been shown (Lenzuni et al. 1991, ApJS, 76, 759) to be the dominant 
positive ion, in a zero metallicity gas at low temperature and intermediate to high 
density. It therefore affects both the number of free electrons and the opacity of 
' the gas. The most recent Hjj~ partition function (Neale and Tennyson, 1995, ApJ, 

454, L169) is an order of magnitude larger at 4000 K than all previous partition 
functions, implying that H^" is a more important electron donor than previously 
O 1 thought. 

Here we present new Rosseland mean opacities for a hydrogen-helium gas of 
1000 < T(K) < 9000 and -14 < log 10 (p [g cm" 3 ]) < -2. In the calculation of 
Oh" these opacities we have made use of the latest collision induced absorption data 

as well as the most recent Hjj~ partition function and line opacity data. It is 
shown that these updated and new sources of opacity give rise to a Rosseland 
mean opacity for a hydrogen-helium gas which is in general greater than that 
calculated in earlier works. 
X ■ The new opacity data are then used to model the evolution of low mass 

0.15-0.8 M Q zero metallicity stars, from pre-Main Sequence collapse to Main 
Sequence turn-off. To investigate the effect of H^ on the evolution of low mass 
zero metallicity stars, we repeat our calculations neglecting as a source of 
electrons and line opacity. We find that H;j" can have an effect on the structure 
and evolution of stars of mass ~0.5 M Q or less. A gray atmosphere is used for 
the calculation, which is sufficient to demonstrate that Hjj~ affects the evolution 
of very low mass stars to a greater degree than previously believed. 
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1. Introduction. 

It has long been speculated that zero metallicity, Population III stars exist, but as yet 
none of these stars have been identified. At the current epoch low mass stars are the most 
numerous stars in the Galaxy, however the first stars are often thought to have formed with 
exclusively high masses. For example Tohline (1980) estimated the minimum mass of a Pop- 
ulation III star to be a massive 1500 M , more recently de Araiijo & Opher (1989) estimated 
the minimum mass to be 50 M . On the other hand Kashlinsky & Rees (1983) proposed 
a method by which stars of mass as low as 0.2 M Q could be formed from primordial mate- 
rial. Recently using hydrodynamical simulations Nakamura & Umemura (2001) suggested 
a bi-modal Population III initial mass function with a minimum mass of between 1 and 
2 M . The discovery of the extremely metal deficient star HE0107-5240 (Christlieb et al. 
2002) which has an estimated mass of 0.8 M Q and [Fe/H]=-5.3 has added credence to the 
possibility that low mass primordial stars formed. All Population III stars with mass less 
than about 0.75 M should still be on the Main Sequence and thus could, if ever observed, 
provide an estimate of the age of the Universe. 

The accuracy of stellar formation and evolution calculations is, among other things, 
dependent upon available opacity data. In the absence of metals the only contributors to 
the opacity of the atmospheres Population III stars are the various hydrogen and helium 
species and ions, and electrons. In a zero metallicity gas at densities common to stellar 
atmospheres it is H~ that dominates opacity at temperatures between about 3500 K and 
7000 K. Below about 3500 K it is collision induced absorption (CIA) by H2-He and H2-H2 
that dominates the opacity. Many authors (Lenzuni et al. 1991; Alexander & Ferguson 1994; 
Saumon et al. 1994; Rohrmann 2001; Allard et al. 2001) have commented upon how H^" can 
act as an electron donor in low metallicity stars, thereby increasing the abundance of H~ and 
hence the opacity. However none have yet investigated the effect of H3 rotation-vibration 
line absorption on a hydrogen-helium gas opacity. 

The zero metallicity opacity calculations of Lenzuni et al. (1991) used the partition 
function of Patch (1968), whereas the opacity calculations of Alexander & Ferguson (1994), 
Saumon et al. (1994) and Bergeron et al. (1995) used the partition functions of Chandra et 
al. (1991). The partition functions of Patch (1968) and Chandra et al. (1991) are an order 
of magnitude smaller at T~4000 K, than the latest most accurate H3 partition function 
of Neale & Tennyson (1995). The Neale & Tennyson (1995) partition function was 
computed by direct summation of 200 000 rotation- vibration energy levels. This most recent 

partition function has been used extensively in the modeling of cool stellar, white and 
brown dwarf atmospheres (Bergeron et al. 1997; Hauschildt et al. 1999; Saumon & Jacobson 
1999; Rohrmann 2001; Allard et al. 2001). The cool stellar and brown dwarf atmospheres of 
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Hauschildt et al. (1999) have subsequently been incorporated in evolution models of brown 
dwarfs and low mass stars by Chabrier and Baraffe (1997) and Baraffe et al. (1997, 1998, 
2002). In particular Baraffe et al. (1997, 1998) studied very low mass, low metallicity stars 
of [M/H]> —2.0 and Saumon et al. (1994) studied zero metallicity very low mass stars with 
the Chandra et al. (1991) partition function. However neither the Hjj~ partition function 
of Neale & Tennyson (1995) or any Hg linelist has been used in evolution calculations of 
low mass, zero or very low metallicity stars. Recent very low and zero metallicity low 
mass stellar evolution calculations by VandenBerg et al. (2000), Marigo et al. (2001), and 
Siess et al. (2002) use the low temperature (T< 10 000 K) Alexander & Ferguson (1994) 
opacities. Whereas Richard et al. (2002) use low temperature opacities based on Kurucz 
data, see Turcotte et al. (1998) and Profhtt (1994). The minimum mass to which these very 
low metallicity evolution calculations extend down to is 0.5 M Q . As we show below, zero 
metallicity stars of mass less than about 0.5 M are affected by Hg acting as an electron 
donor. 

In this paper we present new low temperature Rosseland mean opacities for a zero 
metallicity gas in LTE. These opacities have been calculated both with and without the 
Neale & Tennyson (1995) partition function and the Neale et al. (1996) linelist. We 
then incorporate these new opacity data within the CESAM stellar evolution code (Morel 
1997) and examine the effect of H3 on the pre-main and Main Sequence evolution of zero 
metallicity stars. 

2. Low temperature opacity. 

For compatibility with the CESAM stellar evolution code (Morel 1997), the Rosseland 
mean opacity must be determined from the state variables temperature (T), mass density 
(p) and hydrogen mass fraction (X). The opacity of a gas is dependent upon its chemical 
composition, so the abundances of all the species within the gas or plasma must first be 
determined. For this purpose thermodynamic and chemical equilibrium are assumed, which 
for low densities allows the uses of the Saha equation to determine the number densities 
of individual species. By combining these number densities with the cross sections of the 
relevant scattering or absorption process we can obtain the total opacity for the process in 
the gas. Rosseland mean opacities are then computed as below. 
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2.1. The equation of state. 

The species that are considered in this work are listed in table 1 together with the source 
of the rotational- vibrational partition function if relevant. We use an iterative approach to 
solve for chemical equilibrium using the Saha equations for each chemical species. 

For the initial conditions of our calculation a neutral H2 and He gas is used, the number 
densities of all other species are set to zero. The initial number densities of H2 and He are 
determined from temperature, p and X. For most of the values of T and p considered in 
this work H 2 , H or He are the most abundant species. So we first determine the dominant 
hydrogen species by calculating a first order number density of H from H 2 using the Saha 
equation: 

n(H 2 ) Ck<fr)rx<m (1) 



n(H) 2 Qt(H) 2 ^ \ kT 
where n(x) is the number density of the species x, Qt{%) is the total partition function of 
species x, and x(H 2 ) is the dissociation energy of H 2 . The total partition function, Qt(x), 
which for H 2 is the product of translational, rotational-vibrational and electronic partition 
functions. In this work we neglect electronic states higher than the ground state and set the 
electronic partition function equal to the degeneracy of the ground electronic state. 

Equation 1 is solved iteratively using a Newton-Raphson technique and forcing the 
conservation of H nuclei. Using the Saha ionization equation we then determine the first 
order number densities of H + , H~, He + , H^ , H^, HeH + and e~. Again an iterative Newton- 
Raphson technique is employed, taking care to conserve charge and hydrogen and helium 
nuclei. 

The determination of the number density of H3 is slightly more difficult. The usual 
formation reaction of H^ in the interstellar medium is: 

H+ + H 2 ^ H+ + H (2) 

see for example Herbst. (2000). In chemical equilibrium the number density of H^ is often 
many orders of magnitude smaller than that of H^ which in turn is many orders of mag- 
nitude smaller than the number densities of H 2 and H. This can result in problems with 
machine accuracy. Instead we estimate the number density of H3 from H 2 or H, depending 
upon whether H 2 or H is more abundant. For instance if H 2 is dominant then we use the 
hypothetical reaction: 2H 2 ^ H3 + H + e~. This reaction has the Saha equation: 

"( H *) 2 = [Qt(H 2 )] 2 ,x\ (3) 

n(H+)n(H)n(e") Q T (R+)Q T (R)Q T (e~) P \kT ) 1 ] 

where x is the difference in energy between products and reactants, this is equal to the 
dissociation energy of H 2 (4.52 eV) plus the ionization potential of H (13.60 eV) minus the 
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proton affinity of H 2 (4.38 eV), so 13.74 eV. This approach increases the stability of 
the calculation of the number density of hht . It is valid because, by definition, in chemical 
equilibrium each species is in equilibrium with every other species. 

Once the first estimate for the number density of H^ has been obtained, the process is 
reiterated, using the new number densities, until converged values for the number density of 
each species are obtained. 

Figure 1 shows an illustration of the temperature dependences of the number densities 
of each species at p = 1CT 6 g cm" 3 and X=0.72. Importantly at p = 1CT 6 g cm -3 Fhj" is 
the dominant positive ion below 3500 K, so will have a strong effect on the abundances 
of ions and electrons below this temperature. Figure 2 shows an illustration of the mass 
density dependences of the number densities of each species at T=3500 and X=0.72. Here 
Hjj" becomes the dominant ion above p ~ 2 x 10~ 6 . Figure 3 shows the number density 
of H and H^ calculated using the Hjj" rotational-vibrational partition functions of Neale 
& Tennyson (1995), Chandra et al. (1991) and Patch (1968) as a function of temperature 
with p = 10~ 6 and X=0.72. Clearly the number densities of H3 and H are strongly 
dependent upon which partition function is used, indeed above 4000 K use of the Neale & 
Tennyson (1995) partition function yields an order of magnitude more H3 than the use of 
other partition functions. However at this mass density H + becomes the dominant ion above 
3700 K and the effect of H^ on H~ quickly diminishes with increasing temperature. As H~ 
is an important source of opacity, the opacity is sensitive to which partition function is used 
in the calculation. 

The limits of the partition functions used in this work effectively determine the tem- 
perature limits for which this equation of state is reliable. For all but H^ the partition 
functions extend up to at least 9000 K, and down to 1000 K. For H^ the Neale & Tennyson 
(1995) partition function extends to 8000 K. At the densities of interest to us Hjj" is relatively 
unimportant between 8000 and 9000 K. Therefore our equation of state is reliable between 
1000 to 9000 K. The main short coming of the Saha equation is that it does not account for 
pressure ionization which occurs at very high densities. This imposes an upper density limit 
above which our equation of state is not valid. Mihalas et al. (1988) show that for T less 
than 10 000 K, the upper density for which the Saha equation is valid is around logp = —2, 
we take this as the upper density limit for our equation of state. 
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2.1.1. The effect of metals on the number density of H£ . 

Hjj" is readily destroyed in a dissociative recombination reaction with electrons. Hence 
high abundances of electrons act to reduce the abundance of H^. Many species of metals 
have an ionization energy that is significantly lower than that of hydrogen, so that, in a gas 
containing metals, it is the metals that govern the abundance of free electrons. In this way 
metals have a strong effect on the abundance of H^. As all known stars contain some metals 
we briefly investigate the effect that these metals have upon Hj. 

The equation of state is dependent upon the relative abundance of each metal. As the 
stars in which H3 is an important electron donor are likely to be metal poor we use the metal 
abundance mix found in the very metal poor star HE0107-5240, discovered by Christlieb et 
al. (2002). Where the abundances of these metals are unmeasured for HE0107-5240 we use 
the metal abundances based upon a theoretical Population III supernova yield calculated by 
Umeda & Nomoto (2003). In our equation of state we include what are likely to be the ten 
most abundant metals in HE0107-5240: C, N, O, Ne, Na, Mg, Si, S, Ca and Fe. We included 
the first ionization states, and the anions of the metals within our equation of state. 

Figure 4 shows the number densities of Hjj~, H + , e~ and the positive metal ions (Z + ) as 
function of temperature, at a density of 10 -6 g cm -3 . From this figure, Hjj - is the dominant 
positive ion in a gas between temperatures of 3000 and 3600 K with a elemental abundances 
similar to HE0107-5240. Figure 5 shows the number densities of Hj|", H + , e~ and the positive 
metal ions (Z + ), as function of metal mass fraction. We have used the same metal mix as 
for figure 4 with a constant hydrogen mass fraction of X=0.72 and a density of 10~ 6 g cm" 3 . 
The total number density of the positive metal ions only exceeds the number density of 
H3 for metal mass fractions above 10~ 3 . It should be noted that the metal mix found in 
HE0107-5240 has a very high carbon to iron ratio relative to solar ([C/Fe]= 4.0), but has 
a comparatively low abundance of more readily ionized elements such as Na ([Na/Fe]= 0.8) 
and Fe. Hence this metal mix will result in fewer free electrons at low temperatures than a 
solar metal mix. For a solar metal mix (Grevesse & Sauval 1998) with the same temperature 
and density, we find that the metals are the main electron donors for metal mass fractions 
above about 10~ 6 . The relative importance of H^ as an electron donor is not only dependent 
upon metallicity, but also the relative abundances of the metals within the mix. 

Density also has a role in determining the dominant positive ion, high densities favor 
Hjj" low densities favor the metal ions. Hence H^ is more likely to have an effect on very 
metal poor low mass dwarf stars rather than low mass giants. 
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2.2. The opacity. 

The frequency dependent continuous opacity can now be calculated using the number 
densities of the various species. In this work the set of subroutines developed by Booth & 
Lynas-Gray (2002) were used to calculate the continuous opacity contributions from H I, 
H~, He - , He I, He II bound-free and free- free, H^ free-free, Rayleigh scattering of He I, H 2 
and H I and Thompson scattering by e~. We have also included the Lebedev et al. (2000) 
H^ free-free and bound-free opacity data. The method used by Booth & Lynas-Gray (2002) 
for each absorption or scattering process is summarized in table 2 along with all the other 
opacity sources. We have not included any photodissociation processes for HeH + as these 
are considered unimportant or important only at extreme ultra violet wavelengths (Roberge 
and Dalgarno 1982). We have also included H I line absorption by using the STARK 
subroutine from the ATLAS12 (Kurucz 1993) program. H I line absorption is important for 
temperatures greater than ~7500 K. 

In this work we make use of the most recent CIA data, which are in a tabular form. 
We use the H2-He CIA data of J0rgensen et al. (2000), which covers the temperature range 
1000 to 7000 K and the wavenumber range 25 to 20 000 cm -1 . The H2-H2 data of Borysow 
et al. (2001) which covers the temperature range 1000 to 7000 K and the wavenumber range 
20 to 20 000 cm -1 . The H-He data of Gustafsson and Frommhold (2001) which covers the 
temperature range 1500 to 10 000 K and the wavenumber range 50 to 11 000 cur 1 . No CIA 
data for H-He collisions were available to Lenzuni et al. (1991) and Alexander & Ferguson 
(1994), but these data have been used by Rohrmann et al. (2002) in the calculation of 
synthetic white dwarf spectra. Our calculations indicate that at low hydrogen mass fractions 
and logp ~ —8 and logT ~ 3.5 H-He CIA contributes significantly to the Rosseland mean 
opacity and is the dominant source of continuous opacity between wavelengths of ~3.3 to 
~13 /im. We have neglected H-H 2 CIA, however, since it has only a small effect on opacity in 
a narrow temperature range when H 2 and H are of near equal abundances. At low hydrogen 
mass fractions, low temperatures and moderate to high densities He-He CIA may have an 
effect on opacity. We have not found any He-He opacity data, so have not been able to 
include this source of opacity in our calculations. 

For the purposes of our calculation the CIA opacity tables of J0rgensen et al. (2000); 
Borysow et al. (2001); Gustafsson and Frommhold (2001) are interpolated using a bicubic 
spline. Frequencies at which CIA is an important source of opacities are covered by the 
frequency ranges of the data tables. Beyond the upper temperature limit of the CIA data 
tables many opacity sources exceed, by several orders of magnitude, the opacity from CIA. 
However CIA from H 2 -He and H 2 -H 2 collisions remains important for some temperatures 
below 1000 K. So the calculation of the total opacity is limited only by the lower temperature 
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limits of the H 2 -He and H 2 -H 2 CIA data tables, which is 1000 K. 

To investigate the effect of direct absorption by H3 the ab initio rotational-vibrational 
linelist of Neale et al. (1996) was used. This linelist contains over 3xl0 6 lines between and 
15 000 cm -1 all of which were used. A thermal Doppler line profile was used to distribute 
the intensity of each line, this approximation to the true line profile is valid as all the H3 
lines are weak such that the extended line wings contribute little to the overall opacity. 
was found to contribute up to 15% of the Rosseland mean opacity via line absorption. This 
is far smaller than the indirect effect has on opacity through electron donation, which 
can increase the Rosseland mean opacity by over a factor of 3. However it implies that 
lines may well be observable in cool Population III stellar spectra. We have identified a 
mistake in figures 1 and 2 of Neale et al. (1996). An incorrect degeneracy factor of 2 instead 
of 8/3 was used to calculate the absorption by hotter lines in these figures, which resulted in 
the quoted absorption and emission being reduced by up to 25%. The mistake affects only 
figures 1 and 2 and has no bearing on the Einstein A coeflicents of the Neale et al. (1996) 
linelist. 

The ion HeH + has a strong electric dipole moment and so will have a strong infrared 
spectrum. We have not been able to find an extensive set of HeH + line intensities, so we 
cannot include HeH + line absorption into our opacity calculation. The line opacity of 
only contributes significantly to the Rosseland mean opacity when the opacity of the gas is 
low and the number density of is high. The number densities of HeH + exceed only 
at high temperatures, when the Rosseland mean opacity is high. We would expect the total 
line absorption intensities of HeH + per molecule to be weaker than that of . So we believe 
that HeH + will have little effect on the opacity, but this needs to be verified. 

Figure 6 shows the frequency dependent opacity of the dominant contributors to the 
total continuous opacity at T=3500 K, p = 10~ 6 g cm -3 and X=0.72. The total opacity is 
calculated by summing all the contributors to the opacity. This figure illustrates the strength 
of CIA, H~ and H3 line absorption, hence the importance that these sources of opacity 
are correctly accounted for. The H3 line absorption has a clear effect on the frequency 
dependent opacity between 1500 and 8000 cur 1 . This indicates that line absorption may 
well be observable in certain cool dense metal free stars, brown and possibly cool white 
dwarfs. 

The Rosseland mean opacity is defined by: 
Where B u is the Planck function, k r is the Rosseland mean opacity and k v is the monochro- 
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matic opacity. The denominator is evaluated analytically, but the numerator must be evalu- 
ated numerically. We use 30 000 integration points, which is sufficient to account for H3 line 
absorption and to converge the integration. The partition functions we have used to calcu- 
late the number densities of the various species limit us to temperatures above 1000 K. We 
therefore only give opacity data for temperatures above 1000 K. This minimum temperature 
is sufficient for our stellar models. 

Figure 7 show the Rosseland mean opacity calculated in this work along with those 
of Lenzuni et al. (1991) and Alexander & Ferguson (1994) at p = 10~ 7 and X=0.72. The 
Alexander & Ferguson (1994) data were interpolated using a bi-cubic spline. A density of 
p = 10~ 7 is too low for the effect of on the number density of H~ to be evident. However 
at this density H3 line opacity is strong enough to contribute 13 % of the Rosseland mean 
opacity at 3250 K. The Rosseland mean opacity of this work for T< 3500 K is significantly 
higher than those of Lenzuni et al. (1991) and Alexander & Ferguson (1994). At these 
temperatures the opacity is dominated by CIA and a little H3 line opacity. The older works 
of Lenzuni et al. (1991) and Alexander & Ferguson (1994) do not include line opacity 
and use older CIA data which is weaker than the more recent CIA data, see Borysow et 
al. (2001) and J0rgensen et al. (2000). So a greater opacity at these temperatures is to be 
expected. 

Figure 8 shows a plot of the Roseland mean opacity at p = 10~ 5 g cm -3 and X=0.72. 
The opacities were calculated with the partition functions of Patch (1968), Neale & 
Tennyson (1995), Chandra et al. (1991) and also without accounting for FLj". The Rosseland 
mean opacity of Lenzuni et al. (1991) is also plotted, but that of Alexander & Ferguson (1994) 
does not extend to these high densities. Again the difference between the opacity of Lenzuni 
et al. (1991) and this work are evident, at these temperatures collision induced absorption 
is important for T< 4500 K and accounts for part of the difference. H3 line absorption is 
weak making up ~1.1% of the Rosseland mean opacity at 4000 K. But the most important 
factor is clearly attributable to the effect of on the opacity. The partition functions of 
Patch (1968) and Chandra et al. (1991) lead to an under prediction of the abundance of 
H and hence the opacity. It would appear that has an effect on the opacity several 
times stronger, and to about 500 K hotter, than predicted by the partition functions of 
Patch (1968) and Chandra et al. (1991). At higher densities still, these differences are even 
larger and electron donation becomes even more important. The contribution of line 
absorption to the Rosseland mean opacity peaks at about logp = —7.5 and T=3125 K, at 
which it is 15 % of the opacity. 

We give Rosseland mean opacities for a range of densities, temperatures and hydrogen 
mass fractions in table 3. 
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3. Stellar evolution calculations. 

We use the ID CESAM stellar evolution code (Morel 1997) for our calculations. The 
CESAM code has been widely applied to Helio- and asteroseismology, but has also been used 
for studies low mass stars (0.6 M Q ) (Morel and Baglin 1997). 

The basic input physics is identical to that described by (Morel 1997), except for the 
nuclear reaction rates, opacities and equations of state, which are as follows. We use the 
recent NACRE (Angulo et al. 1999) compilation of nuclear reaction rates, which had been 
previously incorporated into the CESAM code and used for Solar evolution models (Morel et 
al. 1999). The recently improved and extended OPAL equation of state (Rogers & Nayfonov 
2002) which we refer to as OPAL02 is used for all temperatures and densities, from the 
core to the surface of the star. The most recent OPAL opacities Iglesias & Rogers (1996) 
which we refer to as OPAL96, are used for temperatures above 9000 K and the opacity data 
discussed in the previous section is used for T < 9000 K. The equation of state discussed 
in section 2.1 is used solely to calculate the number densities for use in the calculation of 
the low temperature opacities (T < 9000 K). The OPAL02 equation of state provides the 
various thermodynamic quantities necessary for the calculation. To investigate the effects of 
H3 on these models we calculate models using 2 different low temperature opacity functions. 
Firstly we perform a set of calculations using the Neale & Tennyson (1995) H^ partition 
function and secondly we repeat the calculation while completely neglecting hLt . 

It should be noted that the that the OPAL96 opacities are limited to values of log R < 
1, where R = p/(T 3 x 10~ 18 ). As the mass of a star decreases its density increases, such that 
for low metallicity stars of about 0.7 M and below, values of logi? exceed 1 and can exceed 
2 in a star of 0.4 M . In such circumstances we have extrapolated the OPAL96 opacity 
tables to higher values of logi?. An extrapolation to logi? = 2 requires an extrapolation of 
an order of magnitude in density or an extrapolation of just over a factor of 2 in temperature. 
However the region in which we are forced to extrapolate the OPAL96 opacities corresponds 
to the convective regions of the stars. As convective, and not radiative, energy transport 
is dominant the extrapolation of the OPAL96 opacities is unimportant. To test that this 
extrapolation is unimportant we perturbed the extrapolated OPAL96 opacity by multiplying 
by logi? when logi? >1. We found that the calculations performed with and without this 
opacity perturbation were identical. 

The minimum mass of star that we study is imposed by the boundaries of the OPAL02 
equation of state tables. These equation of state tables do not cover all the temperatures 
and densities within our model stars of mass less than about 0.15 M . Hence we limit our 
calculation to masses of 0.15 M and above. 
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A gray atmosphere is used for the surface boundary condition, the temperature strati- 
fication follows a T(t) relationship based on Hopf's law. A thorough discussion of this and 
other model atmospheres used in CESAM can be found in Morel et al. (1994). In the zero 
metallicity model atmospheres of Saumon et al. (1994), significant differences between the 
gray and non-gray (log(g) = 5) model atmospheres do occur for T e fr < 4500 K. Consequently 
for stars below this T c g- gray atmospheres give relatively poor boundary conditions. Non- 
gray atmospheres effectively increase the average atmospheric opacity relative to the gray 
case (Saumon et al. 1994). In this work the inadequacies of the gray atmosphere are rela- 
tively unimportant, because this is primarily a differential study in which we seek to show 
that the use of the H3 partition function of Neale & Tennyson (1995) causes to have a 
greater effect on stellar evolution than previously thought. To enable us to calculate more 
realistic evolutionary models, we are in the process of calculating non-gray atmospheres for 
zero metallicity very low mass stars. 

The initial conditions for our calculation are those of a collapsing pre-Main Sequence 
star at the top of the Hayashi track, as described in Morel (1997). Briefly the pre-Main 
Sequence star is in quasi-static equilibrium and energy is released only through gravitational 
collapse. The star is fully convective so that entropy is constant throughout the star. The 
initial conditions are described by a contraction constant, c. For a 1 M star Morel (1997) 
suggests a contraction constant of 0.02 L©Mq 1 K _1 . In this work we use a constant value of 
0.015 LqM^K- 1 , we find that the value of the contraction constant has only a small effect 
the position of the resultant star upon the Main Sequence. 

Our models are computed without mass loss and rotation. We follow Marigo et al. 
(2001) and use an initial hydrogen mass fraction of X=0.77. We set the metal mass fraction 
to be Z= 10 -14 , this is low enough to make the effect of metals negligible, whist avoiding 
numerical problems with extremely low and zero values of Z. This leaves a Helium mass 
fraction (Y) of marginally under 0.23. 

4. Results and Discussion. 

Figure 9 shows the evolution paths of zero metallicity model stars between 0.4 and 0.15 
M calculated with and without the effects of H3 . Symbols are placed on the curves when 
the core hydrogen mass fractions are 10~ 2 , 10~ 3 and 10~ 4 . Initially the pre-Main Sequence 
stars of 0.4 M Q and less, with and without H3 contract along identical paths, the density 
within the atmosphere at this stage is too low for to have any effect. As the star becomes 
more dense and before the stars reach the Main Sequences the paths separate as starts 
to become the dominant positive ion within the atmosphere. At zero age main sequence the 
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star with Hg has a lower T e g- and luminosity, than the star in which Hg is not accounted 
for. This trend continues throughout the Main Sequence lifetime of the stars. For stars 
with higher masses of 0.5 and 0.45 M Q we find that affects the late pre-Main Sequence 
and early Main Sequence evolution, but that the evolutionary tracks on the Main Sequence 
converge as the star becomes hotter in the late main sequence, and H^ is destroyed. Our 
evolution tracks are in agreement with the zero metallicity calculations of Marigo et al. 
(2001) which cover masses greater than 0.7 M Q and also with Siess et al. (2002) which cover 
masses greater than 0.8 M . 

This tendency for model stars of less than 0.5 M with to have lower T e ff and 
luminosity than their counterparts without H^ is shown again in figure 10. This figure 
shows the 14 GYr isochrone for model stars calculated with and without the effects of Hg , 
with symbols placed at intervals of 0.05 M . The three curves represent models calculated 
using the Neale & Tennyson (1995) Hg partition function, the Chandra et al. (1991) Hg 
partition function and also neglecting ELt . For 14 GYr old stars of 0.15 M , the difference 
in T e g between models calculated with the Neale & Tennyson (1995) partition function and 
those in which is neglected is over 200 K and the difference in luminosity is about 18 
%. The models calculated with the Neale & Tennyson (1995) partition function show that 
H^ has a significantly greater effect, which extends to higher masses, than that predicted 
by models calculated with the Chandra et al. (1991) partition function. This illustrates the 
importance of using the Neale & Tennyson (1995) partition function in low mass evolution 
models. The luminosity and effective temperature of the 0.2 M and 0.15 M zero metallicity 
Main Sequence model stars of Saumon et al. (1994) are also plotted in figure 10. As Saumon 
et al. (1994) used the Chandra et al. (1991) partition function, their results would have 
been expected to lie upon the isochrone calculated with the Chandra et al. (1991) partition 
function. However the models of Saumon et al. (1994) lie closer to the isochrone calculated 
with the Neale & Tennyson (1995) partition function. Saumon et al. (1994) used non-gray 
atmospheres whereas in this work we have used gray atmospheres, so the reason for this 
difference is likely to be due to the break down of the gray approximation in very low mass 
stars. Saumon et al. (1994) commented that using non-gray models effectively increases the 
atmospheric opacity over gray model atmospheres. Hence non-gray stars have a lower T c g- 
than their gray counterparts. 

It is clear from figure 10 that the strength of the effect of Hg increases with decreasing 
mass. As long as H remains the dominant source of opacity the strength of the effect of 
Hg will continue to increase with decreasing mass. The effect of Hg within the atmosphere 
will diminish when collision induced absorption becomes the dominant source of opacity, at 
about T e fj < 3000 K. However in a shell below the photospheres of such stars hotter denser 
gas exists in which H~ remains the dominant source of opacity and Hg will again affect 
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opacity. 

Figure 11 shows the evolution of the central temperature for model stars of mass 0.8, 
0.6, 0.4, 0.25 and 0.15 M , calculated with and with out the effects of Hg . Again it is seen 
that Hg has no effect on the the central temperatures of model stars of mass 0.6 M and 
above. For the model stars of 0.25 and 0.15 M the increased atmospheric opacity due to 

acts to reduce the central temperature. This in turn reduces the nuclear reaction rates 
and so increases the hydrogen burning lifetime of the model star. also acts to reduce 
the pre-Main Sequence and early Main Sequence central temperatures of the 0.4 M model. 
However for most of the Main Sequence the central temperature of the 0.4 M model star 
which neglects Hjj - is cooler than the 0.4 M model which accounts for H^. This results in 
a longer hydrogen burning lifetime for the model which neglects Hjj". The reason for this 
is that the 0.4 M model star without Hg is fully convective at zero age Main Sequence. 
Whereas the 0.4 M model star which accounts for H^ starts zero age Main Sequence with a 
radiative interior, but the model star which neglects H^ starts zero age main sequences with 
a fully convective interior. Convective mixing in a fully convective star makes the hydrogen 
throughout the star available for fusion in the core. Conversely in a star with a radiative 
interior there is only very limited mixing of material, so that only the hydrogen in the core 
is available for fusion. Hence the convective model star without H^ has a longer lifetime 
than the radiative model star with Hg , even though the former is more luminous. It appears 
that Hg within the atmosphere, can have an effect upon the structure and energy transport 
mechanisms within these low mass stars that have interior temperature gradients close to the 
adiabatic temperature gradient. So slight changes in the temperature gradient brought about 
by increased atmospheric opacity can change the main energy transport method. All the 
stars of mass less than 0.4 M have a fully convective interior at zero age main sequence, all 
of the stars develop a radiative interior at some point later in their main sequence evolution. 
After developing a radiative interior the rate of evolution of the star along the main sequence 
increases dramatically. 



5. Conclusion. 

We have calculated new low temperature zero metallicity Rosseland mean opacities 
using the latest collision induced absorption opacity data and with the latest H^ partition 
function and line opacity data. The results of this calculation confirm that at high densities 
and intermediate temperatures H^ becomes the dominant positive ion. Furthermore using 
the Hg Neale & Tennyson (1995) partition function instead of an older partition function, 
causes the temperature and density ranges at which Hg forms to expand as well as increasing 
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its abundance by up to an order of magnitude. This in turn increases the abundance of 
electrons and hence H~ and other negative ions, thereby affecting the opacity of the gas by 
a greater degree than was previously thought. The line opacity at certain temperatures 
and densities can contribute up to 15 % of the Rosseland mean opacity. H^" lines may well 
be observable within cool Population III stars and all models of such stars should account 
for direct absorption by IPj". 

The new Rosseland mean opacities have been included in evolution models of very low 
mass zero metallicity stars, enabling us to investigate the effect of Hjj" on the evolution of 
these stars. For late pre-Main Sequence and early Main Sequence stars of mass < 0.55 M Q 
we find that Hjj~ can cause the atmosphere to cool and the star to drop slightly in luminosity. 
For stars of 0.45 M Q and below this effect becomes strong enough to affect the Main Sequence 
structure and lifetime of the star. 

The results presented here show that as the mass of a zero metallicity star is reduced 
the effect of Hjj" upon its evolution, structure, luminosity and temperature increase. As the 
Neale & Tennyson (1995) partition function is much larger than the partition function of 
Chandra et al. (1991), used by Saumon et al. (1994), its corresponding effect on opacity is 
larger. Hence if the effect of H3 continues to grow in objects of below 0.15 M , then Hjj" will 
effect the evolution of zero metallicity very low mass stars and brown dwarfs too a greater 
extent than previously believed. 

Our anonymous referee is thanked for their constructive criticism which has helped to 
improve this paper. We thank Dr. P. Morel for use of his stellar evolution code CESAM 
and Dr. A. Borysow for making the results of her CIA calculations publicly available. We 
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Fig. 1. — The number densities of the species considered in the calculation as a function of 
temperature, at p = 1CT 6 g cm' 3 . The solid curves represent negatively charged species, 
the dotted curves represent positively charged hydrogen species, the dashed curves represent 
neutral hydrogen species and the dot-dashed curves represent helium containing species. 
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Fig. 2. — The number densities of the species considered in the calculation as a function 
of density, at T=3500 K. The solid curves represent negatively charged species, the dotted 
curves represent positively charged hydrogen species, the dashed curves represent neutral 
hydrogen species and the dot-dashed curves represent helium containing species. 
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Fig. 3. — The number densities of Hjj" and H as a function of temperature, at p =1CT 6 g 
cm -3 . The three curves were calculated with the Hjj" partition functions of Neale & Tennyson 
(1995) solid curves, Chandra et al. (1991) dashed curves and Patch (1968) dotted curves. 
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Fig. 9. — The evolutionary tracks of 0.4, 0.3, 0.2 and 0.15 solar mass zero metallicity stars 
calculated by both neglecting and including the effect of H3 upon opacity. Symbols are 
placed at values of the hydrogen mass fraction of 10~ 2 , 10~ 3 and 10 -4 at the center of the 
star. 
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Table 1: The species used in the calculation, the source of their rotational- vibrational parti- 



tion funct 


ions is quoted where applicable. 


Species 


T range (K) 


Reference 


H 2 


1000-9000 


Sauval & Tatum (1984) 


H 






H 






H+ 






H 2 + 


1500-18000 


Stancil (1994) 


H 2 - 


1000-9000 


Sauval k Tatum (1984) 


H 3 + 


500-8000 


Neale & Tennyson (1995) 


He 






He+ 






HeH+ 


1000-9000 


Sauval & Tatum (1984) 


e~ 
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Table 2: The sources of continuous monochromatic opacity used in the calculation, with 
temperature and wavenumber limits if applicable. 

Limits 
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He bound-free 








Kurucz (1970) Method 


T-Tp fvpp.fypp 








Knrnc7 f!970 N ) Method 


He + bound-free 








Kurucz (1970) Method 


He + free-free 








Kurucz (1970) Method 


H bound-free 








Method of Gray (1976) 


H free-free 








Method of Gray (1976) 
using Gaunt factors from 
Karzas & Latter (1961) 


free-free 


658 - 


28 500 


> 1400 - 10 000 


Bell (1980) 


H+ b-f, f-f 


4500 - 


- 10 000 


2000 - 46 000 


Lebedev et al. (2000) 


H 2 Rayleigh scat. 






> 9.76 


Kurucz (1970) method. 


H Rayleigh scat. 






< 82 000 


Kurucz (1970) method. 


He Rayleigh scat. 






> 17.2 


Kurucz (1970) method. 


c~ Thomson scat. 








Method of Gray (1976) 


H3 Line opacity 






< 15 000 


Neale et al. (1996) 


H I Line opacity 








Subroutine from ATLAS 12 
Kurucz (1993) 
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Table 3: Rosseland mean opacities (kr) for a range of hydrogen mass fractions (0 < X < 1), 
densities (lCT 14 < log 10 (p [g cm" 3 ]) < 10~ 2 ) and temperatures (1000 < T(K) < 9000 ). The 
complete version of this table is in the electronic edition of the Journal. The printed edition 
contains only a sample. 



X 



logi (p fe cm 3 D T ( K ) logi («fl [cm 2 g x ]) 



0.70E+00 
0.70E+00 
0.70E+00 
0.70E+00 
0.70E+00 
0.70E+00 
0.70E+00 
0.70E+00 
0.70E+00 
0.70E+00 
0.70E+00 
0.70E+00 
0.70E+00 
0.70E+00 
0.70E+00 
0.70E+00 
0.70E+00 
0.70E+00 
0.70E+00 
0.70E+00 
0.70E+00 



-5.50E- 
-5.50E- 
-5.50E- 
-5.50E- 



-00 
-00 
-00 
-00 



-5.50E+00 
-5.50E+00 
-5.50E+00 
-5.50E+00 
-5.50E+00 
-5.50E+00 
-5.50E+00 
-5.50E+00 
-5.50E+00 
-5.50E+00 
-5.50E+00 
-5.50E+00 
-5.50E+00 
-5.50E+00 
-5.50E+00 
-5.50E+00 
-5.50E+00 



2.00E+03 
2.25E+03 
2.50E+03 
2.75E+03 
3.00E+03 
3.25E+03 
3.50E+03 
3.75E+03 
4.00E+03 
4.25E+03 
4.50E+03 
4.75E+03 
5.00E+03 
5.25E+03 
5.50E+03 
5.75E+03 
6.00E+03 
6.25E+03 
6.50E+03 
6.75E+03 
7.00E+03 



-4.158E+00 

-4.181E+00 

-4.189E+00 

-4.179E+00 

-4.095E+00 

-3.802E+00 

-3.374E+00 

-2.973E+00 

-2.593E+00 

-2.157E+00 

-1.738E+00 

-1.362E+00 

-1.031E+00 

-7.399E-01 

-4.811E-01 

-2.486E-01 

-3.780E-02 

1.552E-01 

3.335E-01 

4.994E-01 

6.551E-01 



